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We investigate direct above-threshold ionization in diatomic molecules, with particular emphasis 
on how quantum interference is altered by a driving field of non-vanishing ellipticity. This inter- 
ference may be either temporal, i.e., related to ionization events occurring at different times, or 
spatial, i.e., related to the electron emission at different centers in the molecule. Employing the 
strong-field approximation and saddle-point methods, we find that, in general, for non-vanishing 
ellipticity, there will be a blurring of the temporal and spatial interference patterns. The former 
blurring is caused by the electron velocity component perpendicular to the major polarization axis, 
while spatial interference is washed out as a consequence either of s — p mixing, or of the temporal 
dependence of the ionization prefactor. Both types of interference are analyzed in detail in terms of 
electron trajectories, and specific conditions for which sharp fringes occur are provided. 

PACS numbers: 33.80.Rv, 33.80.Wz, 42.50.Hz 



I. INTRODUCTION 

Atoms and molecules interacting with a strong laser 
field may be ionized by absorbing more photons than 
necessary. This highly nonlinear phenomenon is called 
above-threshold ionization (ATI) and has attracted in- 
creasing attention since its 1979 discovery by Agostini 
and co-workers [l|. Recently, impressive progress has 
been achieved in the study of ATL For example, great 
emphasis has been placed on employing ATI as a tool for 
the attosecond imaging of molecular orbitals [2,l3l- Imag- 
ing applications are based on the physical mechanism be- 
hind this phenomenon, in which the active electron (i) is 
released in the continuum, (ii) is accelerated by the field, 
and (iii) reaches the detector without returning or after 
rescattering with its parent molecule. The former and the 
latter scenario are known as direct ATI, and as rescat- 
tered or high-order ATI (HATI), respectively. Quantum 
mechanically, this electron may be released along several 
pathways, and the corresponding transition amplitudes 
will interfere. 

Previous studies of ATI have shown that for the strong- 
field ionization of molecules, there are two kinds of quan- 
tum interference displayed in the ATI spectrum. In the 
first scenario, for a given energy, the electronic wave 
packet may reach the continuum at events occurring at 
different times in one cycle of laser field [4] . Hence, the 
transition amplitudes associated with ionization along 
different orbits in the time domain will interfere. This 
leads to a fringe structure in the ATI spectrum [3, Q, 
which is present regardless of whether the target is an 
atom or a molecule. For that reason, we will refer to this 
type of interference as "temporal interference" . Recently, 
ATI spectra with unusual two-photon separation in the 
direction perpendicular to the laser polarization [7| have 
been revealed due to the temporal interference. 

In the second scenario, due to the multi-core structure 
of molecules, there may be electron emission from spa- 



tially separated centers, which results in a multi-slit like 
interference pattern in the ATI spectrum. In the con- 
text of linearly polarized fields and diatomic molecules, 
spatial-interference effects have been widely studied [3- 
Il6j . Such studies range from simplified molecular mod- 
els to exact three-dimensional time-dependent calcula- 
tions. Employing the molecular strong-field approxima- 
tion (SFA), a generalized spatial interference condition 
for diatomic molecules that accounts for the orbital ge- 
ometry and s — p mixing has been derived [14l4lq . For 
such fields, however, the dynamics are mainly confined 
along the field polarization axis. 

Elliptically polarized laser fields exhibit a nonvanish- 
ing electric-field component perpendicular to the major 
polarization axis, which affects the motion of the elec- 
tron in the laser field. Hence, they have been proposed 
as a resource for controlling strong-field phenomena by 
varying the ellipticity. For example, it is known that el- 
liptical polarization has great significance as a tool to 
produce isolated attosecond pulses ^T^ and to increase 
the contribution of long orbits in the HATI process jl^l ■ 
So far, the ellipticity dependence of ATI has been widely 
studied since the mid 1990s. Most studies, however, fo- 
cus on the influence of the field ellipticity on the angular 
dependence of the ionization rate of ATI |19l - l22| . while 
the imprints left on the ATI spectrum by a nonvanishing 
driving-field ellipticity have received little attention. 

In an early paper Q, features related to the above- 
mentioned temporal interference have been identified ex- 
perimentally for direct ATI. Therein, it was found that 
the photoelectron yield as a function of the driving-field 
ellipticity behaves in two very distinct ways for HATI 
and direct ATI energy ranges. While in the former case 
this yield decreases monotonically for larger driving field 
ellipticity, following the predictions of classical models, 
for direct ATI the electron ellipticity distributions dis- 
play several "shoulder" structures. These structures have 
been related to the quantum temporal interference be- 



tween the two sets of orbits that contribute to ATI in 
this energy region. However, how the elhpticity influ- 
ences the two orbits stiU remains unexplored to a great 
extent. 

In this work, we are interested in how the driving-field 
ellipticity influences the above-mentioned two kinds of 
interferences in the ATI spectra of diatomic molecules, 
exemplified by N2 and Ar2. Throughout, we focus on di- 
rect ATI, for which temporal interference effects are also 
prominent. Employing the molecular SFA and saddle- 
point methods, we find that, in general, the temporal in- 
terference in the ATI spectra will become blurred in the 
elliptically polarized field. For electron momenta along 
the major polarization axis, however, the temporal inter- 
ference remains clear and the positions of the interference 
minima shift towards lower photoelectron energies with 
increasing ellipticity. These features are well explained in 
terms of the ellipticity dependence of the initial electron 
velocity along each orbit, and we show that the downshift 
of the interference minima is responsible for the "shoul- 
der" structure observed in the ellipticity distributions of 
ATI electrons in Ref. ^. 

Furthermore, for the spatial interference, clear inter- 
ference minima which do not depend on the ellipticity 
are found for N2 if only s or p states are taken. How- 
ever, if s — p mixing is introduced, sharp fringes are only 
present for linear polarization if the molecule is aligned 
perpendicular to the field axis. For all other cases, these 
minima are blurred even for linear polarization. For Ar2, 
the ATI spectra display several clear minima caused by 
the spatial interference. The reasons for the differences 
between N2 and Ar2 molecules are discussed in detail in 
terms of the shapes of their highest occupied molecular 
orbitals (HOMO). In this view, we propose a necessary 
condition to observe the spatial interference in the ATI 
spectrum in experiments. 

This paper is organized as follows. In Sec. II, we pro- 
vide a brief discussion of the model, including how to ob- 
tain the ATI amplitude and two-center interference con- 
ditions. Subsequently, we present the ATI spectra of the 
N2 and Ar2 molecules in the elliptically polarized laser 
field, and discuss how quantum interference is altered by 
the driving field of non-vanishing ellipticity. Finally, in 
Sec. IV our conclusions are given. Atomic units (a.u.) 
are used throughout unless otherwise indicated. 



where 



1 1"°° 



(2) 



is the semiclassical action. In Eqs. ([TJ and ([2]), Ip gives 
the ionization potential, and E(t) and A{t) denote the 
laser electric field and the vector potential, respectively. 
For simplicity, we assume that all the molecular structure 
is contained in the ionization prefactor 



(p + A(t)|-r.E(t)|V^o> 



(3) 



This is the most widespread assumption when performing 
the SFA modelin g of strong-field phenomena in molecular 
systems (|13l. [2a - |29{ : see also the reviews [30,[31|), and 
implies that the active electron leaves from the geomet- 
rical center of the molecule. In this work, we employ the 
single-active orbital approximation, i.e., we assume that 
the initial state |'0o) is the HOMO. This orbital is written 
as the linear combination of atomic orbitals (LCAO), and 
the core is assumed to be frozen. This latter approxima- 
tion is reasonable for molecules with heavy nuclei [32| . 
Explicitly, 



V'o(r) = X1' 



V'a(r+5:) + (-l) 



Ic-r 



"V'aCr-— ) 



(4) 
where the wavefunctions V'a(r) give the atomic orbitals, 
R is the internuclear distance, l^ is the orbital quantum 
number, rua is the magnetic quantum number and A^ = 
ma for gerade (g) and Aq = m+l for ungerade (u) orbital 
symmetry [33J. 



Saddle-point equation and temporal 
interference 



For sufficiently high intensity and low frequency of the 
laser field, the temporal integrations in the amplitude 
^ can be mathematically evaluated by the saddle-point 
method with high accuracy [41, l24|, i.e., we seek solutions 
such that the action ^ is stationary. For the specific 
case of an elliptically polarized field confined to the xy 
plane, the saddle-point equation to be solved reads as 



[Px+A^{t)]'^ [Py + Ay{t)] 



/p = 0. 



(5) 



II. MODEL 

A. Transition Amplitude 

Within the molecular SFA, the transition amplitude 
for direct ATI from the initial bound state \tjJo) to a final 
Volkov state with drift momentum p is given by |4l.l24ll25| 

/oo 
dt{p + A{t) \-r-E{t)\ Vo)exp[z^d(p,i)], 
-00 

(1) 



where Px and Py denote the final momentum in the x 
direction and in the y direction, respectively. Physically, 
Eq. ^ ensures the conservation of energy at the ioniza- 
tion time t. In terms of the solutions tg of Eq. ([S]), the 
transition amplitude ([T} can be written as 

Md{p) (x^Cs{ts)dionexp[iSdip,ts)], (6) 

S 

where the prefactors 



CsiQ 



2m 



d^Sd{p,u)/dti 



(7) 



are expected to vary much more slowly than the action at 
each saddle for the above-stated approximation to hold 
[35l |. For each final momentum p, there are two solutions 
of Eq. ([5|) within one period of the laser field, whose 
contributions will interfere. This interference is present 
regardless of whether the target has one or more cen- 
ters, i.e., for atoms and molecules, and has been exten- 
sively studied in the literature for linearly polarized driv- 
ing fields [1143 • 

In this work, we employ the monochromatic elliptically 
polarized laser field 



E(t) = 






with vector potential 



A{t) 



v/T+? 



■.{xsinut ~ ^ycosut), (8) 



:{xcosu}t + ^ysinujt), (9) 



where ^ is the laser ellipticity. In this case, all cycles are 
equivalent and the amplitude Mdip) can be intuitively 
understood as the coherent superposition of the contri- 
butions of the two orbits, which may add constructively 
or destructively. For vanishing electron momenta, these 
orbits will start at two consecutive maxima of the driving 
field E(i). These times will approach each other as the 
momentum increases. Throughout, we will refer to the 
orbits starting in the first and second half cycle of the 
field dH) as orbit 1 and orbit 2, respectively. For a discus- 
sion of these solutions for linearly polarized monochro- 
matic fields see [H-Isi El. 



C. Spatial interference condition 

Due to the two-center structure of the diatomic 
molecule, one may expect electron emission from spa- 
tially separated centers in the molecule. This results in 
minima and maxima in the ATI spectra, which, in the 
present model, are implicit in the ionization dipole matrix 
element dion- Using the length gauge with field-dressed 
molecular bound states [30|, the ionization prefactor ^ 
can be rewritten as 



ia-rjj.Q+Ao 



dion = -^Ca[exp{ip ■ R/2) + {-1) 



exp(-ip ■ R/2)]E(t) • idpMp) 



(10) 



where tpaip) is the Fourier transform of the atomic or- 
bitals V'a(r) given by Eq. ^ with p = p-l-A(i). The two- 
center interference is related to the terms in the bracket 
in Eq. pil|) . This generalized interference condition ac- 
counts for s — p mixing and the orbital geometry, and 
reduces to simpler interference conditions if only s or p 
states are present. For example, the terms in bracket 
reduce to 2 cos(p • R/2) for s states and gerade (g) sym- 
metry. Thus, interference minima are present if 



where k denotes an integer number. Similarly, the inter- 
ference minima are given by the condition 



p ■ R = 2fc7r 



(12) 



for p states and gerade symmetry. If, on the other hand, 
the HOMO has ungerade (u) symmetry, the above-stated 
scenario is reversed, i.e., the destructive interference con- 
dition for p states is given by Eq. PT|) . while Eq. p^ 
holds for s states. 

In this paper, the molecular axis is kept in the plane 
xy spanned by the driving-field ellipticity. The frame 
of reference (x ,y ,z ) of the molecule is rotated by an 
angle rj in this plane, defined with regard to the major 
polarization axis of the laser field. This means that x = 
sin rjy + cos rji , y — cos jyy — sin ryx , and z = z . 

The wave functions "00(1") in Eq. Q will be approxi- 
mated by a Gaussian basis set 

Mr) = E ^-(^')^^' (y')^"' (^')''^' ^M-Cy), (13) 

V 

where the coefficients Cq, by^ I3x',y',z' and the expo- 
nents Ct, are extracted from the quantum chemistry code 
GAMESS-UK [H]. Throughout, we employ a 6-31G ba- 
sis so that only s and p states are considered. For a 
Gaussian basis set, the wave function ?Aq(p) in momen- 
tum space can be written as 



V'a(p) =^Ktpy{p). 



(14) 



Usually, the explicit formula of (pv{p) is complicated. 
However, if only s and p states are used in the orbital 
construction, ipv{p) is simplified as 

3/2 

<^.(P) = i-VxY" I "^3/2+;^ exp[-pV(4a)], (15) 

where /„ = 0, 1 and the subscript x indicates the orien- 
tation of the p states employed to construct the orbital. 
For instance, a a orbital would require p states along the 
X axis, i.e., x — x' y whereas a n orbital would require p 
states along the y or z axis, i.e., x — v' ot z' . 

In this work, we will study the N2 and Ar2 molecules, 
both of which have a HOMOs. In combination with Eqs. 
(dni and ([5]), the prefactor in Eq. ^TU\\ can be specified 
as 

„ _ T> 

dion = i C0S{——) E Caby[pxEx{t) +PyEy{t)]'Ds + 



a.(s),v 



,p R 



sin(— — ) E Caby{[{pl-2(y) cos ri+pxPySmT]]x 



a.(p).v 

Ex{t) + [{pl - 2Cy) smrj+pxpy cos r]]Ey{t)}Vp 



(16) 



p-R= (2fc + l)7r, 



(11) 



where P, = e^Tr^/VC^^^ and Vp = e"^ 77^/^(2(1^^). 
Note that, in Eq. (1161) . the momentum components relate 
to the electric field frame of reference. 
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FIG. 1: (Color online) ATI spectra of the N2 molecule, coded 
in false colors in the (Ep, 9) plane and plotted in a logarithmic 
scale. In panels (a) and (d), both orbits 1 and 2 were included, 
for 5 = and ^ — 0.36, respectively, while in panels (b), (c), 
(e) and (f) the contributions of individual orbits were taken. 
In panels (b) and (c), the field has been chosen to be linearly 
polarized (C = 0), while in panels (e) and (f) the ellipticity is 
5 — 0.36. The major axis of the polarization ellipse is along 
the molecular axis, i.e., for alignment angle r; = 0. The laser 
intensity and wavelength are 7 = 2 x 10^* W/cm^ and A =800 
nm, respectively. 



III. RESULTS AND DISCUSSION 

In the results that follow, we will analyze how the 
above-mentioned temporal and spatial interference are 
influenced by a non-vanishing driving field ellipticity. 
This analysis will be performed in terms of electron or- 
bits. For the temporal interference, we will focus on the 
N2 molecule, which has been widely studied in the lit- 
erature, both theoretically and experimentally. For all 
parameters employed in Sec. IIIIAl the structure caused 
by electron emission at separate centers gets washed out, 
so that temporal interference may be clearly assessed. In 
our studies of spatial interference effects (Sec. IIIIB[) . we 
consider the N2 and Ar2 molecules. Such molecules have 
very different HOMO geometries and equilibrium inter- 
nuclear distances. 



A. Temporal interference 

Figure [T] shows the ATI spectra of the N2 molecule, 
coded in false colors in the {Ep, 9) plane, where 9 — 
Z(R, p) denotes the angle between the internuclear dis- 
tance and the drift momentum, and Ep — p^/2 is the 




FIG. 2: (Color online) Time evolution of the two components 



of the vector potential A(f) 
and the electric field E(t) 



= (xcosa;f + ^ysinojt) 



= (x sin Ldt — ^y cos tut) , with 



5 = 0.36 [panels(a) and (b), respectively]. The two points 
ti and t2 denote the times associated to the two ionization 
events 1 and 2 that occur within a field cycle, for an electron 
with the same final momentum p. These events are indicated 
by the colored circles in the figure. For simplicity, we have 
defined the time related to a field crossing Exit) = for the 
electric-field component along the major polarization axis as 
i = 0. This time also gives a maximum of the larger vector 
potential component Ax(t). Note that Ax{t\) = Ax{t2) as 
both events are equidistant in time from i = 0. 



kinetic photoelectron energy upon reaching the detector, 
in units of the ponderomotive potential Up. We assume 
that the major axis of the polarization ellipse is along the 
molecular axis, i.e., that the alignment angle rj = ^(R, x) 
is vanishing. For linear polarization [see Fig.[IJa)], many 
clear minima can be observed in the spectra. These min- 
ima are due to the destructive interference between orbit 
1 and orbit 2. To illustrate this more clearly, we calculate 
the ATI transition probability density |Mrf(p)p employ- 
ing only individual orbits, for the same parameters as in 
Fig. [Ha). These results are displayed in Figs. [TJb) and 
(c), and unambiguously show no minima. 

With increasing ellipticity, the interference minima be- 
come blurred and, for most angles 9, completely disap- 
pear as the ellipticity reaches £, — 0.36, as shown in Fig. 
[Tfd). In fact, only the minima along and in the vicin- 
ity of the major polarization axis {9 — and tt) remain 
visible for this ellipticity. Sharp minima are only present 
exactly along the axis |43| . 

This behavior is due to the fact that, in general, a non- 
vanishing driving-field ellipticity affects the initial elec- 
tron velocity along orbit 1 and orbit 2 unequally. The 
initial velocity of the electron at the ionization time t 
is p -t- A(i), and thus the velocity along the minor po- 
larization axis is given by Vy = pshi9 + Ay{t). In Fig. 
Ufa), we present the time evolution of the vector poten- 
tial A(t) with ellipticity ^ = 0.36. The two times ti and 
^2 represent the ionization times related to orbit 1 and 



2, respectively, and are equidistant from the electric field 
crossing Ex{t) — [see Fig. [H^b)]. The figure clearly 
shows that the directions of the vector potential Ay are 
opposite at the two ionization times. Thus if the electron 
leaves along orbit 1, the absolute value of its initial veloc- 
ity Vy should be larger than if it leaves along orbit 2. This 
holds except for linear polarization, i.e., if Ay(t) = 0, or 
for electron momenta along the major polarization axis, 
i.e., for 6 — or TT. The ionization rate at a specific ion- 
ization time t is evaluated by w(i, wq) ~ WQ{t)wi{t,vo) 
[ioll. where 



Wo{t) = 



Wi{t,Vo) 



4(2/p)^ 
E{t) 

E{t)TT 



exp 



exp 



3E{t) 



i^Ip) 



3/2 



E{t) 



mp) 



1/2 



(17) 



(18) 
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FIG. 3: (Color online) ATI spectra of the N2 molecule along 
the major polarization axis as a function of the alignment 
angle rj, coded in false colors in the {Ep, -q) plane. EUipticity 
^ = for (a) and £, — 0.36 for (b). The parameters of the 
laser field are the same as in Fig. \T\ 



and vo 



Vy is the initial velocity. Hence, the 



larger the initial velocity vq is, the lower the ionization 
rate at that ionization time will be. This implies that, 
for elliptical polarization, ionization along orbit 2 will be- 
come larger than ionization along orbit 1. Consequently, 
the destructive interference between two electron orbits 
will disappear in the ATI spectra. Here we assume that 
the transverse velocity of the electron is Py — psind < 0. 
li Py >Q the above-stated results are reversed. 

This is in agreement with Figs.[lje) and (f), in which 
we present the individual contributions to the ATI spec- 
tra from orbit 1 and orbit 2, respectively, for ellipticity 
^ = 0.36. In general, the probability density |Mci(p)p as- 
sociated with orbit 1 is much lower than that related to 
orbit 2. This is in stark contrast to what is observed for 
linear polarization, for which the two probability den- 
sities are identical [see Figs, [ijb) and (c)]. However, 
for drift momentum p along the major polarization axis 
{6 = and vr), the ATI probability densities are the same, 
because the two orbits exhibit the same initial perpen- 
dicular velocity \vy\ = \Ay{t)\. Therefore, the coherent 
superposition of the corresponding transition amplitudes 
will lead to clear interference minima in the ATI spectra, 
as observed in Fig. [T{d). 

Next we will further study the interference minima in 
the ATI spectra along the major polarization axis. In 
Fig. 131 we calculate the ATI spectra of the N2 molecule 
with = as a function of the alignment angle rj. The 
most eye-catching features in the spectra are the clear 
minima due to the destructive interference between the 
two electron orbits. With increasing ellipticity, the ATI 
spectra seem to be squeezed and the positions of the min- 
ima are moving down to lower energies. For example, in 
Fig.[2Ia) the position of one of the minima is about 2.hUp 
for ry = and ^ = 0, but it becomes 1.9t/p in Fig. ^h) 
for ^ = 0.36. Furthermore, Fig. [3] shows that the down- 
shift of the minima with the field ellipticity is much more 
dramatic for higher photoelectron energies than in the 
low-energy region. 

In order to understand the behavior of the minima in 
Fig. [21 we calculate the phase difference between the two 
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FIG. 4: (Color online) (a) Phase difference between the con- 
tributions of the two orbits in Eq. ((6| for parallel align- 
ment r; = and momentum along the major polarization 
axis {e = 0) with C = (dashed lines) and C = 0.36 (full 
lines); action (black), prefactor (red), and total transition 
amplitudes (green). The gray arrows indicate the values of 
the energy when the total phase difference is equal to Stt and 
Stt for ^ = and ^ = 0.36. (b) Saddle points associated with 
the solutions of Eq. ([5]) for the two orbits as functions of the 
ellipticity from ^ = to ^ = 0.36, for parallel alignment, mo- 
menta along the major polarization axis, and kinetic energies 
Ep = Q.^Up (black full lines), Ep = Up (red dashed lines) and 
Ep = 2Up (blue dotted lines). The gray arrows indicate the 
real values of the saddle points for ^ = and f = 0.36. 



orbits in Eq. ^ , which determines the positions of the 
interference minima in the ATI spectra. In Fig. JH^a), 
we exhibit the phase difference between the amplitudes 
M^^^(p) and M^^^(p), the actions S'd(p, ti) and S'd(p, ^2), 

and the prefactors Ci(p, ti)^^-^^ and C2(p,i2)'^ion associ- 
ated with the two orbits for parallel alignment ry = 0, 



and field polarizations ^ = and ^ = 0.36 (dashed and 
full lines, respectively). The figure shows that the phase 
difference between Sd{p,ti) and Sd{p,t2) (black lines) 
is smaller for ^ — 0.36 than for ^ — 0. This is caused 
by the real parts of two saddles ii and i2 approaching 
each other for elliptical polarization [43], as illustrated 
in Fig. |4l(b). Moreover, the approaching of the two sad- 
dles becomes more significant in the high-energy than in 
the low-energy region of the photoelectron spectra [see 
the gray arrows in Fig. IDJb)]. Therefore, the decrease of 
the phase difference associated with the action is much 
more dramatic for higher photoelectron energies than in 
the low-energy region. For the prcfactor, the phase dif- 
ference changes much less markedly with regard to the 
ellipticity, as shown in the red lines in Fig. lU^a). This is 
in agreement with the assumption of the steepest descent 
method that the prefactor varies more slowly with time 
than the action [3a]. Thus, the change in the total phase 
difference [green lines in Fig. |4l(a)] with the ellipticity is 
mainly determined by the action. When the difference in 
the total phase is equal to (2fc -|- l)7r, fc = 0, 1, 2, ..., de- 
structive interference between the two orbits will occur. 
In Fig. mja), the positions of the total phase difference 
with fc = 1 and fc = 2 are marked with arrows. The two 
positions are about 2.5Up and 1.3Up for linear polariza- 
tion and move down to about 1.9Up and l.lt/p, respec- 
tively, for elliptical polarization. Moreover, the downshift 
of the minima in the high-energy region is more signifi- 
cant than that in the low-energy region. These results are 
consistent with the behaviors of the interference minima 
observed in Fig. |31 

Furthermore, the behavior of the real parts of the sad- 
dle points in Fig. IHJb) can be explained in terms of the 
ellipticity dependence of the initial ionization velocity vq. 
In the classical limit, the real parts Re(ti) and Re(i2) are 
associated with the classical ionization times ii and ^2 of 
an electron in the field ([5]), obtained if /p — > 0. In a lin- 
early polarized laser field, if an electron is released with 
drift momentum p along the major polarization axis, it 
will reach the continuum with the vanishing initial ve- 
locity Wo = 0. In an elliptically polarized field, however, 
the electron should have a nonvanishing initial velocity 
in order to compensate the motion induced by the small 
component of the field. In the classical limit, the total 
initial velocity vo(i) of the electron is perpendicular to 
the laser electric field E(i) at the corresponding ioniza- 
tion time t [13, im, i.e., E(t) • Vo(i) = 0. The condition 
Vo(t) _L E(i), t = ti,t2 requires that Vx{ti) > and 
Vx(t2) > for the parameters employed in this paper. 
Since v^ = p + Aocosujii/y/l + ^^ and Ax{t) > for 
both ti and ^2 [see Fig. [2Ia)], an increase in the value 
of Vx with the driving-field ellipticity requires that ti in- 
creases in order to approach the time for which A^it) 
is maximal, i.e., t — [see Fig. [5]. Moreover, for the 
photoelectrons in the high energy region, the ionization 
time is closer to t = v^. Hence, a larger initial ve- 
locity Vx = —Ey/ExVy — £,^Ao cosijjt/ yjl + ^2 will be 
required to satisfy the condition vq -L E(t). This results 
in a larger increase of the ionization time. The same 
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FIG. 5: (Color online) ATI spectra of the N2 molecule com- 
puted using only the s states [panels (a) and (c)] and the p 
states [panels(b) and (d)] from its HOMO, coded in false col- 
ors in the {Ep, 9) plane. The major axis of the polarization 
ellipse is along the molecular axis, r] = Q. Ellipticity ^ = 
for (a) and (c) and ^ = 0.36 for (b) and (d). In all pan- 
els, the two-center interference minima are indicated by the 
white curves in the figure. The remaining field and molecular 
parameters are the same as in Fig. [T] 



argument applies to orbit 2, with the difference that, be- 
cause this orbit is located after the electric-field crossing 
[see Fig. [UJb)], ^2 should decrease. As a direct conse- 
quence, for increasing ellipticity the two ionization times 
should approach each other, and the approaching of the 
two times becomes more significant for higher photoelec- 
tron energies than in the low-energy region. These results 
are in good agreement with the two saddle-point solu- 
tions as functions of the drivin g-fie ld ellipticity shown in 
Fig.HKb) (or see Fig. 4 in Ref. (l^). 

It is worth mentioning that, because the interference 
minima in Fig. [3] move down to lower energy with in- 
creasing ellipticity, for fixed photoelectron energy there 
will be fluctuations in the ATI spectrum. This interest- 
ing phenomenon has been observed in experiments [231 ] 
and was called the "shoulder" structure in the ellipticity 
distribution of the ATI electrons. In our work, the under- 
lying mechanism of the "shoulder" structure is explained 
in terms of the ellipticity dependence of the initial ion- 
ization velocity. 



B. Two-center interference 

Finally, we will study the spatial interference due to 
electron emission at the two spatially separated centers in 
N2 and Ar2. The former molecule has a 3a'g HOMO with 
a high degree oi s — p mixing. For the latter molecule, 
the HOMO is a 5cr„ orbital, in which the p states domi- 
nate. Furthermore, Ar2 exhibits a very large equilibrium 
internuclear distance {R = 7.2 a.u.), so that, according 
to Eq. pI7)) . several two center minima are expected to 



be present in the momentum range of interest. 

In Fig. [5l we display the ATI spectra of N2, if only 
the contributions from the s or p states to the HOMO 
are taken. These spectra are coded in false colors in the 
(-Ep, 9) plane. The major axis of the polarization el- 
lipse is chosen to be parallel to the molecular axis, i.e., 
at the alignment angle 77 = 0. Fig. EJa) shows that for 
the s-state contributions and linear polarization, apart 
from the minima caused by the temporal interference, 
there are other two minima, whose positions satisfy the 
two-center interference condition, Eq. (|lip , shown by the 
white lines in the figure. With increasing ellipticity [see 
Fig-El^)]: the minima from the temporal interference are 
blurred as discussed in the previous section, but the min- 
ima from the two-center interference are still very clear. 
That is because the two-center interference condition. 
Eg. (|TT]) . is given for field-dressed initial bound states 
[39J. In this case, it is independent on the ellipticity. 
Similarly, for the p-state contributions, the two-center in- 
terference minima appear in the middle of the each panel 
in Figs.[5jb) and (d), because the corresponding electron 
momentum is perpendicular to the molecular axis, i.e., 
p • R = 0. Hence, the interference condition, Eq. ([T^. 
always holds. 
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FIG. 6: (Color online) (a) and (b) Logarithm of ATI spectra 
of the N2 molecule, coded in false colors in the {Ep, 9) plane. 
The major axis of the polarization ellipse is perpendicular 
to the molecular axis, rj = 7r/2. Ellipticity ^ = for (a) 
and ^ = 0.36 for (b). (c) and (d) ATI spectra obtained by 
employing only individual orbits in the (-Bp, 6) for ^ = 0. The 
white curves indicate the minima caused by the two-center 
interference. 



However, when we consider the full 3ag HOMO of the 
N2 molecule and include s — p mixing, the minima from 
the two-center interference can not be observed in the 
ATI spectra. This holds for linearly and elliptically po- 
larized fields, as shown in Figs, [ija) and (d), and even 
if only the contributions of individual orbits are taken, 
as illustrated in Figs. [TJc), (b), (e), and (f). There are 
two main reasons for it. First, ii s—p mixing is included, 
the interference condition also depends on the sum of 
dpipa{p) in Eq- p^ - where p = p + A(i). This influence 
is expected to be considerable for N2, as its Scr^ HOMO 
has a high degree oi s — p mixing. Second, the ioniza- 
tion time t has a positive imaginary part, which can be 
related to the width of the barrier through which the elec- 
tron must initially tunnel Q- Thus, the sum in Eq. pi?)) 
is usually not purely real or imaginary, and therefore, the 
absolute value of Eq. p^ will never completely vanish, 
i.e., no clear minima will be observed in the ATI spectra. 

Nevertheless, there is a specific case in which clear spa- 
tial interference fringes are present in the full ATI spectra 
of N2, ii s — p mixing is included. Fig. ^a.) shows the 
spectra for linear polarization with the alignment angle 
77 = 7r/2. Apart from the minima from the temporal in- 
terference, there are two other minima. In order to show 
them more clearly, we calculate the ATI transition prob- 
ability density |Mrf(p)p by employing only individual or- 
bits to remove the pattern associated with the temporal 
interference. Figs. Efc) and (d) display the ATI spec- 
tra without the temporal interference, where the clear 
minima from the spatial interference are observed. The 
reason for the appearance of the spatial interference in 
this specific case is that for the perpendicular alignment 
r] = t: 12 and for linear polarization, the ionization pref- 



actor in Eq. (fT6|) can be simplified as 

„ . r> 

djon =iC0s(— — ) ^ CahvPxEx{t)Vs + 

ism{——) ^ CcbvPxPyEx{t)Vp, 

Q(p),f 



(19) 



where the momentum py = py in the y direction is real. 
Thus, ipxExit) can be factorized in both the first and 
the second terms in Eq. (|19|) . and the remaining terms 
are time-independent and purely real. Using Eq. (|19l) . 
we obtain the minima from the destructive spatial inter- 
ference in the ATI spectra, shown by the white curves in 
Figs, mja), (c), and (d), which agree well with what is 
observed in the spectra. On the other hand, for elliptical 
polarization, the imaginary parts in the two terms in Eq. 
(dni) are different. Hence, they cannot be factorized and 
the absolute value of Eq. PH)) will never vanish. This is 
also the case if the field is linearly polarized, but rj ^ 7r/2, 
as seen from Eq. (|19p . Fig. Eljb) shows that as the ellip- 
ticity increases to i^ = 0.36, the minima from both the 
spatial and temporal interference in the spectra fade. 

In contrast, the ATI spectra of the Ar2 molecule 
present clear minima from the two-center interference, 
as seen in Fig. [7fa). In addition. Fig. Eljb) shows that, 
with increasing ellipticity, most minima from the tempo- 
ral interference are blurred. However, the spatial minima 
remain unaltered. The critical reason for the difference 
between Ar2 and N2 molecules is that, in the 5cr„ HOMO 
of Ar2, the contributions of the p states dominate and 
the contributions of the s states are practically vanishing 
for the internuclear separation in this work. Hence, the 
corresponding two-center interference condition is given 
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FIG. 7: (Color online) ATI spectra of the Ar2 molecule, with 
internuclear distance R = 7.2 a.u., coded in false colors in 
the (-Ep, 6) plane. The major axis of the polarization ellipse 
is along the molecular axis, j; = 0. The results in panels 
(a) and (b) have been computed for ellipticity ^ = and 
^ = 0.36, respectively. The minima caused by the two-center 
interference are indicated by white curves. The remaining 
field parameters are the same as in Fig. [T] 



by Eq. PT|) . which holds for pure p states and ungerade 
symmetry. This condition, shown by the white lines in 
Fig. [71 is independent of the driving field and is consistent 
with the interference minima in the simulation. Further- 
more, since Ar2 has a very large equilibrium internuclear 
distance, R — 7.2 a.u., several two-center minima are ob- 
served in the momentum range of interest according to 
Eq. dn]). 



IV. CONCLUSIONS 

In summary, we theoretically study temporal and spa- 
tial interference in the direct ATI spectra of aligned di- 
atomic molecules in an elliptically polarized laser field. 
Temporal interference relates to ionization events occur- 
ring at different times within a field cycle, and spatial 
interference is associated with electron ionization at spa- 
tially separated centers in the molecule. Throughout, we 
have employed the molecular SFA, saddle-point methods 
and field-dressed bound states. We have also assumed 
that only the HOMO is active. 

With increasing ellipticity, there will be a blurring of 
the temporal-interference patterns, except when the pho- 
toelectron momenta p are parallel to the major polar- 
ization axis. This blurring is caused by the fact that a 
nonvanishing driving field ellipticity affects the electron 
velocities unequally for the two ionization events occur- 
ring within each field cycle. For linear polarization, these 
events are equally probable. In order, however, to com- 
pensate the motion induced by the small component of 
the elliptically polarized field, it is necessary that, for one 
of the orbits, the electron be released in the continuum 
with a larger velocity. This will lead to a suppression 



of the corresponding ionization probability, and hence to 
a disappearance of the temporal interference patterns. 
Nevertheless, if the electron is released along the major 
polarization axis, the absolute value of the initial velocity 
will be the same for the two ionization events. Hence, the 
temporal-interference minima in the ATI spectrum will 
remain clear. For these specific parameters, we have ana- 
lyzed in detail how a non- vanishing ellipticity shifts these 
minima towards lower energies, and have traced them to 
a decrease in the phase difference between the transition 
amplitudes related to both orbits. These results are in 
agreement with the studies performed in |23| . 

Our results also show that, in general, the spatial, two 
center interference is much more robust with regard to 
the electric field parameters than the temporal interfer- 
ence, as long as the HOMO does not exhibit a high degree 
of s—p mixing. As a consequence of s—p mixing, however, 
or of the temporal dependence of the ionization prefactor 
Eq. (jlOp . the interference patterns become blurred, even 
for linear polarization. Only for the very specific case 
of linear polarization and a molecule aligned perpendic- 
ular to the field direction will sharp fringes be observed. 
This is exemplified by computations in N2 . On the other 
hand, if either s or p states are dominant in the HOMO, 
the patterns will be very clear over a large parameter 
range and practically independent of the driving-field el- 
lipticity. This is shown for the ATI spectra of Ar2, which 
display clear spatial-interference patterns due to fact that 
the p-state contributions are dominant. Several minima 
exist due to the large equilibrium internuclear distance. 

Therefore, one can trace all the interference patterns 
observed, or the absence thereof, to the dynamics of the 
electronic wave packet and the velocity with which it 
is released in the continuum. These patterns can be 
controlled by an appropriate choice of the the field and 
molecular parameters. Concrete examples provided in 
this work include the laser-field polarization, the molec- 
ular targets, the angle between the electron momentum 
and the major polarization axis, and the molecular orien- 
tation with regard to the field. This control is related to 
the energy position and to the sharpness of the patterns. 
We also identify a parameter region, for which clear in- 
terference fringes can be easily observed. We expect this 
information to be useful for future experiments. 
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